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ABSTRACT 

The  degree  of  approximation  achievable  by  piecewise  polynomial 
functions  on  certain  regular  grids  in  the  plane  is  shown  to  be  adversely 
affected  by  smoothness  requirements*  -  In  stark  contrast  to  the  univariate 
situation.  For  a  rectangular  grid*  and  for  the  triangular  grid  derived  from  It 
by  adding  all  northeast  diagonals*  the  maximum  degree  of  approximation  (as  the 
grid  size  1/n  goes  to  zero}  to  a  suitably  smooth  function  is  shown  to  be 
0(n  in  case  we  Insist  that  the  approximating  functions  are  in  C^  . 

This  only  holds  as  long  as  P  <  and  P  <  *  respectively*  with  r 

the  total  order  of  the  polynomial  pieces*  In  the  contrary  case*  some  smooth 
functions  are  not  approximable  at  all.  In  the  discussion  of  the  second  mesh*  a 
new  and  promising  kind  of  multivariate  B-spllne  Is  Introduced. 

AMS  (MOS)  Subject  Classification:  41A15*  41A63*  41A25 

Key  words:  multivariate,  splines,  piecewise  polynomial*  smoothness, 
degree  of  approximation*  B-splines 

Work  Unit  No.  3  -  Numerical  Analysis  and  Computer  Science 

1 )  Mathematics  Research  Center*  University  of  Wisconsin*  Madison,  WI  53706 

2)  Department  of  Mathematics  and  Statistics*  University  of  South  Carolina, 
Columbia,  SC  29208 


1 )  Sponsored  by  the  United  States  Array  under  Contract  No.  DAAG29-80-C-0041 

2)  Supported  by  National  Science  Foundation  Grant  No.  MCS81-01661. 


SIGSIPICMCE  AND  EXPLANATION 


^  One  of  the  important  properties  of  univariate  splines  is  that  in  most 
senses  smooth  splines  approximate  just  as  well  as  do  piecewise  polynomials  on 
the  same  mesh.  This  report  shows  this  to  be  untrue  in  the  multivariate 
setting.  In  particular,  it  details  the  cost  in  approximating  power  one  may 
have  to  pay  for  the  luxury  of  a  smooth  piecewise  polynomial  ^proximant.  In  an 
extreme  case,  piecewise  polynomials  of  total  degree  <  r  on  a  rectangular  grid 

with  all  derivatives  of  order  continuous  will  fall  to  approximate  certain 

•c  cV  = 

smooth  functions  at  all  (as  the  grid  goes  to  zero)  unless  p  is  kept  below 
(r-3)/2  . 


During  the  analysis  of  approximation  on  a  certain  regular  triangular 
grid,  a  novel  kind  of  bivariate  B-spline  is  introduced.  This  B-spline,  in 
contrast  to  the  established  multivariate  B-spline  derived  from  a  simplex,  can 
be  made  to  have  all  its  breaklines  in  a  given  regular  grid.  This  makes  it  a 
prime  candidate  for  use  In  the  construction  of  smooth  multivariate  piecewise 
polynomial  approximation,  and  its  properties  will  be  explored  further. 
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1  2 
C»  da  Boot  and  R.  DeVore 

1.  Introdoctton.  One  of  the  important  properties  of  univariate  splines  is  that  in 
most  senses  smooth  splines  approximate  just  as  well  as  piecewise  polynomials.  As  we  shall 
see,  this  is  no  lon9er  the  case  for  multivariate  splines  where  both  the  smoothness  of  the 
spline  and  the  geometry  of  the  partition  can  have  a  limiting  effect  on  the  order  of 
approximation.  This  type  of  limitation  has  already  been  recognized  in  certain  cases  [9], 
such  as  low  order  splines  on  a  rectangular  grid.  For  example,  there  is  no  effective 
approximation  by  C^^^-cublcs  on  rectangular  grids.  The  purpose  of  this  paper  is  to  give  a 
systematic  study  of  this  and  related  questions.  We  restrict  our  inquiry  to  bivariate 
approximation  but  it  will  be  clear  that  our  techniques  extend  to  higher  dimensions. 

Let  n  a  (w^)  be  partition  of  Into  triangles  and/or  rectangles.  Denote  by 

«■  {»^{n)}  the  corresponding  scaled  partition,  with  »j^(n)  ;»  t^/n  ,  all  1  .  Denote  by 

S  )  the  space  of  splines  of  order  r  and  smoothness  p  on  11  ,  i.e.,  s  e 

r,P  n  *  n 

S  .(11)  if  and  only  If  a  S  C*’ (B^)  and,  on  each  t,(n>  ,  s  Is  a  polynomial  of  total 
r,P  n  1  . . 

«• 

order  r  (I.e.  total  degree  <  r  ).  We  are  interested  In  when  S  .(II  )  is  dense  in 

n-1  r,P  n 

Cq(B^)  .  We  study  this  problem  in  detail  for  two  particular  partitions:  Z  :>  ' 

consisting  of  the  squares  fi,i+1 )x(j,J+1 I  ,  and  the  partition  i  which  results 

whan  each  square  of  Z  is  divided  into  two  triangles  by  introducing  the  northeast 
diagonal.  The  technique  developed  for  these  two  cases  can  be  used  for  more  general 
partitions  as  well. 
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In  Snctlon  2,  «m  study  approxlnatlon  on  £  .  This  nodest  example  already  has  the 

n 

salient  features  of  the  general  problem.  We  show  that  effective  If  and  only 

If  P  <  .  Thus,  roughly  speaking,  only  smoothness  up  to  one-half  the  order  of  the 

polynomial  pieces  Is  allowable  In  this  case.  It  turns  out  that  smoothness  also  affects  the 

m 

rate  of  approxlnatlon  In  that  certain  C^-functlons  can  only  be  approximated  to  within 

.  Thus  any  smoothness  condition  reduces  the  achievable  order  of  approxlnatlon. 
This  should  be  compared  iriLth  the  univariate  case  where  the  full  order  of  approxlnatlon  Is 
achievable  regardless  of  smoothness. 

We  study  approxlnatlon  on  A  In  Section  3.  Here,  the  role  of  the  geometry  of  the 

II 

partition  becomes  more  apparent.  The  splines  effective  if  and  only  If 

2r-4 

p  <  — ^  ,  hence  there  is  a  gain  over  the  case  fr<»  roughly  r/2  to  2r/3  .  It  Is 

clear  from  our  techniques  that  this  Is  due  to  the  fact  that  the  partition  A  has  three 

palruflse  Independent  directions,  vis.  (1,0),  (0,1)  and  (1,1),  whereas  £  has  only  two, 

viz.  (1,0)  and  (0,1).  More  generally.  If  the  partition  D  Is  generated  by  m  pairwise 

independent  vectors  (w)iat  we  mean  by  this  la  made  precise  In  Section  4),  then  s^  p^^n^ 

effective  If  and  only  if  p  <  r-1  -  ,  Thus  the  more  directions,  the  higher  the 

allomble  smoothness,  but  of  course  at  the  expense  of  a  more  complicated  partition. 

The  results  just  described  have  two  con^nents.  First  they  say  that  US^  p^^n^ 

dense  In  Cn(B^)  If  P  Is  too  large.  This  rests  on  the  fact  that  S_  .(D  )  will  not 
Si  &  f  p  n 

contain  splines  of  finite  support  when  P  Is  too  large.  Our  approach  for  this  part  of  the 

problem  Is  more  or  less  the  seuae  In  both  cases  K  •  £  and  n  •  A  .  The  second  )ialf  of  the 

analysis  Is  to  show  that  Us  ^(11)  Is  danse  «rtten  P  Is  suitably  restricted.  This 

n  r,P  n 

requires  the  construction  of  appropriate  approximation  methods.  We  develop  different 

methods  for  constructing  such  approximations  In  the  two  cases. 

For  £  It  Is  easy  to  see  that  S_  S  (£  )  Is  effective  when  P  <  — since  S_ 

*  II  TpP  n  2  *1 

then  contains  the  tensor  products  of  univariate  splines  of  order  04-2  and  smoothness  P  . 

«• 

The  more  difficult  problem  Is  to  show  that  S  approximates  any  C^^-functlon  f  to 
within  0(n”*^*^*^)  .  This  is  done  by  approximating  the  derivative  by  splines  of 
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lower  order  and  saoothnesa.  Integrating  and  making  local  corrections  for  the  approximation 
to  £  . 

Our  construction  of  approximants  from  S  s>  ^  completely  different 

flavor,  being  based  on  extensions  of  the  idea  of  multivariate  B-splines.  These  ideas  extend 

readily  to  more  general  partitions.  Recall  that  multivariate  B-splines  are  defined  by 

certain  cross  sectional  volumes  of  simplices.  If  the  simplex  is  replaced  by  a  more  general 

polyhedron,  the  resulting  function  is  still  a  piecewise  polynomial,  but  now  perhaps  with 

fewer  lines  of  discontinuity.  In  particular,  by  a  suitable  choice  of  the  polyhedron,  we  can 

force  the  discontinuities  of  the  resulting  spline  to  be  contained  in  the  mesh  lines  of 

A  .As  mentioned  in  Section  4,  other  choices  for  the  polyhedra  will  handle  other 
n 

partitions  n  of  R  .  We  feel  that  this  viewpoint  for  constructing  multivariate  splines 
may  prove  to  be  very  useful.  For  example,  some  of  the  standard  finite  elements  can  be 
described  this  way  (31. 

Here  are  some  notatlonal  conventions  used.  Cq(A)  stands  for  the  set  of  r  times 
continuously  differentiable  functions  on  the  set  A  with  compact  support  in  A  .  1*1  is 

the  L^-norm  on  R^  ,  and  1*1 (A)  is  the  h^-norm  over  the  set  A  .  Further, 

Ifl  «“  •£•  J«  max  lD^*'®^fl  and  |£|  t“  max  lD^®'^^fl  . 

^  W*  a+@<r  o+S-r 

|_a _J  denotes  the  largest  integer  no  bigger  than  a  ,  i.e.,  the  "floor"  function,  emd 

a**]  denotes  its  companion,  the  "celling"  function,  which  gives  the  smallest  integer  no 

smaller  than  a  .  Further, 

g  <>  (l,2,...,m>  . 

We  also  need  the  difference  operator  a5  defined  on  the  function  g  by  the  rule 

n 

'  (x»yl  ^  I  (-1)^*^  (p}  (D  g(x+ph,y+qh)  . 
p,q»0 

This  is  the  tensor  product  of  the  univariate  r-th  order  difference  operators  in  the 

coordinate  directions.  Explicitly,  A^  •  6^  •  6^  ,  with 

n  infOi  \0,n; 

hi®  '  ^  g(x+pa,y<-pb)  . 

p-0  P 

-  3  - 


2.  AppraodMtion  on  raetangolar  grid*.  Let  £  tm  {o^^}  with  :• 

(Xj^,Xi^^]x[yj,yj^l]  and  x^  :■  y^  i>  1  .  We  are  Interested  In  the  restriction  of  S 

S  _{Z)  to  SOM  coapact  set. 
tfP 

Lsbm  1.  Every  sSSaS^p(£)  can  be  represented  on 

«»  {(x.y)  t  x,y  >  0) 

as  a  linear  cosJalnatlon  of  the  truncated  powers  in 

T  I-  p(I)  «-  {(x-Xj^)^**(y-y^)^‘*  t  i,j,p,q  >  0|  p4q  <r;^>pifj>0}. 

nrooC.  This  is  proved  by  induction,  siidlar  to  the  proof  that  truncated  power 
functions  are  a  basis  for  univariate  splines.  | | | 


Recall  that  denotes  the  partition  t  scaled  by  1/n  .  Let  x^(n)  i-  y^(n)  :•  1/n. 

®n  «■  . 

iheotwB  1.  is  dense  in  if  and  only  if  P  <  . 

ir*3 

Eroof.  If  p  <  ,  then  contains  all  the  truncated  powers 

(Xiy)  I— »  (x-x.  (n))^*\y-y.  (n)  )**'*''  and  hence  the  space  S  _  of  splines  of 

1  ▼  j  T  n 

coordinate  order  P't’2  and  sooothness  P  .  It  is  known  that  Us  (£  )  is  dense  in 

-  n  p+2,P  n 

2 

Cq(K  )  >  see  for  exaaple  (51  . 

Suppose  that  P  >  euid  that  f  8  C  (R*)  with  If  -  s  I  »  o(l)  as  n  ♦  »  ,  for 

«  u  n 

SOM  splines  ^  •  Suppose  without  loss  of  generality  that  f  has  its  support  in 

2  r“3 

R^  .  Since  P  >  — 2~  ,  the  only  truncated  power  functions  in  T  (as  Introduced  in  Lemia  1 ) 

are  of  the  fore  t(x,y)  »  xP<y-yj)^‘*  or  t(x,y)  «  (x-Xj^)^Py**  with  p  +  q  <  r  .  Since 

A|^t  -  0  for  such  t  ,  we  have  from  Lemma  1  that  A^^s  -  0  on  R^^  for  all  s  8  S  .  This 

implies  that  we  also  have  A^s  •  0  on  R.^  for  n-1,2,,..,  from  which  it  follows  that 

n  n  ^ 

A^f  -  0  ,  If  we  divide  by  h^“^  and  take  the  limit  as  h  ♦  0,  we  find  that  «  0  . 

Since  there  are  C^^-functions  for  which  D^'^'’^^f  0  ,  we  have  proved  the  theorem.  ||| 
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Theorem  1  shows  that  only  splines  with  smoothness  less  than  about  one-half  the  order 

ir~3 

r  will  be  effective  for  approximation.  It  turnn  out  that  even  when  p  <  »  the  order  of 

approximation  is  affected  negatively  by  smoothness.  More  precisely,  we  now  show  that  the 

optimal  order  of  approximation  achievable  with  splines  of  order  r  and  smoothness  P  is 

-r+p+1 

n  • 

ir*3 

Theorem  2.  Let  P  <  — ^  and  k  :»  r-p-1  .  Then 

(i)  there  are  functions  f  e  C*(B^)  for  which  dist(f,S|^)  o(n"’‘)  ,  n-**  > 

( ii )  for  each  f  a  Cq<R^)  ,  dlst(f,S^)  »  0{n~'')  . 

Proof.  Me  will  show  that  if  f  8  c“(*^)  and  distCf.S^j)  -  o(n“’'),  then 
vanishes  at  0  ,  and  this  shows  (i).  Actually,  with  a  finer  analysis,  we  could  show  that 
all  k-th  order  derivatives  of  vanish  on  all  of  . 

Suppose  that  f  8  Cq(R^)  and  that  there  are  functions  Sj^  8  such  that 

(2.1)  If  -  8  I  ■»  o(n“’')  as  n  . 

n 

Let  m  be  a  positive  integer  and  set  h  :<«  1/m  .  If  n  =>  bm  with  b  an  integer,  then  u^, 

4^s  is  in  S_  and,  with  g  s=»  A^f  ,  we  have 

n  n  u 

(2.2)  Ig  -  u^l  -  o(n"’')  . 

Now  any  truncated  po%fer  t(x,y)  -  (x  -  or  t(x,y)  <=  xP(y  -  Xj  )+**  in  T 

r  2 

-  T  _(£)  is  annihilated  by  A.  .  Hence,  on  R.  ,  u_  is  a  linear  combination  of  the 
r^p  R  ▼  n 

splines  A^t  ,  with  t(x,y)  -  (x-Xj^(n))^i’(y-yj(n))^'*  in  and  i,j  >  0  .  Such  t  has 

P  <  p,q  and  p4q  <  r  ,  and  therefore  p,q  <  k  .  Thus  u^  is  a  polynomial  of  coordinate 
order  k  on  each  ®j^j^n)  S  with  i,j  >  0  . 

For  given  n  >  0  ,  choose  b  so  that  1/((b+1)m)  <  kn  <  1/(bni)  .  Then  the  points 
(in,0)  ,  0  <  i  <  k  ,  are  in  a^jj(bm)  and  so  “  0  .  Using  (2.1),  we  have 

•^*’(n,0)’^‘°’l  ■  •^*’(n,0)‘^""bm’^'°’l  ^  =  o(T|'')  , 

therefore  D^*‘'®^g(0)  ■  0  .  If  we  now  let  m  ♦  •  and  recall  that  h  «  1/m  and  g  =  A^f  , 
we  find  >f  (0)  -  0  .  as  desired. 


M«  now  prov«  (11).  Without  loss  of  ganornllty  w«  assume  from  now  on  that  f  Is 
supported  In  the  unit  cube  Qi«(0,  IJ^.If  P“  follows  from  the  fact 

that  contains  the  teneor  product  splines  of  order  )c  •  r-p-1  .  For  general  p  ,  the 
argument  Is  more  Involved.  We  need  a  certain  subspace  of  ,  the  space 


those  s  fl  for  which  s 
{x**/**  :  p+q  <  r  ,  p,q  <  k} 


I  .  .  agrees  with  a  polynomial  In  the  span  of 
|o^^(n)  ’ 

,  all  l.j  .  We  will  prove  by  Induction  on  P  the  following 


Claim.  For  any  P  and  any  r  with  P  <  — j-  there  Is  const  so  that,  for  any  f  0 
Cq(»^)  supported  In  Q  and  any  n  ,  there  Is  an  e  0  p*^n*  "uPPOgt*^  the  square 

(P+1)^Q  for  which 


If  -  el  <  const  Ifl  n 
r 


which  In  turn  gives  (11). 

Certainly  the  claim  Is  true  when  p  «  -1  ,  the  case  of  piecewise  polynomial 

approximation  (see,  e.g. ,  [5]).  Suppose  then  that  the  claim  )ms  been  established  for  all 

ir*3  IT  2 

p  <  Pg  and  consider  P  ■  P^  .  Take  r  so  that  P  <  -j-  .  Take  f  0  Cq(R  )  supported 
In  Q  ,  and  let  g  t>  .  Then 

-  /q/q  dCdn  for  all  (x,y)  0  . 

By  Induction  hypothesis,  we  may  choose  u  0  S_  ,  „  ,(S_)  so  that  the  support  of  u  Is  In 
P^Q  and 

(2.4)  Ig  -  ul  <  const  2  "  2)+(p  1  )+1  ^  const  n  • 

Our  approximation  s  to  f  Is  gotten  by  Integrating  u  and  luklng  local  corrections 
using  B-spllnes. 

For  p  -  (p^,P2)  »  IPt  Mp(x,y)  j-  (w)Rp2(y)  r  with  the  univariate  B-spllne 

with  knots  Xj(n),  ...,  Rj  of  order  P+1  and  smoothness  P-1  ,  and 

Is  supported  on  the  Interval  (x^(n),  follows  that  Mp  la  In 

Sr  2  p-i^^n^  *  assume  to  be  normalised  to  )tave  Integral  1. 

Consider  the  spline  Np  j«  R(p4,^)p  •  supported  on  the  square 

«p  ’■  '*(P+l)p/"’'  *(P+l)(p,+1)‘"”  *  ^'(P+DPj'"*'  y(P*1)(P2+1)‘"” 
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and  therefore  the  have  disjoint  support. 

For  each  p  ,  define  a^  (g  "  u)  .  Prom  (2.4)>  it  follows  that 

P  +  1  mV  — 1 

(2.5)  ja^l  <  iQpl  const  •9*p_2  *  const  n  . 

Consider  now  the  spline 

V  :«  u-fCa  N  . 

P  P  P 

For  each  p  ,  we  have 

(2.6)  /q  tg  -  v)  -  0  . 

Hence,  if  we  set 

3(x,y)  :»  v(C,n)  dCdn  , 

then  s  will  be  in  S  )  and  if  (x#y)  6  then 

r,p  n  P 

f(x,y)  -  3(x,y)  -  /Q/^Cg  -  v)(5,n)  dtdn  =  /g  (g  -  v) 

with 

E  :=  (tO,x]x(y^p^^jp^(n),y))  ( tx^p^^  (n),x)x[0,y) )  _ 

Now,  for  any  (x,y)  e  (P+1 )(}  ,  we  )iave  |e|  <  const  n  ^  .  Further,  /  N  is  nonzero  for 

E  p 

at  most  0(n)  values  of  p  .  Therefore,  (2.4),  (2.5)  give 

lf(x,y)  -  a(x,y)l  -  l/g(g  -  u  -  t  ah))  4  »g  -  ul  Ie)  +  max  la  J  Sj  N 

(2.7)  P  .V 

<  const  If^y  n  . 

This  proves  (2.3)  for  the  approximation  just  constructed.  Also,  using  (2.6)  and  the  fact 

2 

that  f  and  u  have  support  in  p  Q  ,  it  follows  that  s  -  0  outside  of 
(p^  +  (p+1))Q  ^  (P+1)^Q  ,  as  desired.  ||| 

3.  Approximation  on  triangnlar  grids.  He  now  show  how  the  results  of  the  last  section 
can  be  extended  to  triangular  grids.  He  focus  on  the  partition  A  which  is  gotten  by 
dividing  each  square  of  Z  into  two  triangles  by  adding  the  northeast  diagonal.  It  will  be 
clear,  however,  that  our  techniques  apply  to  more  general  partitions,  a  point  made  in 
Section  4.  He  begin  by  developing  methods  for  constructing  smooth  spline  approxlmants.  This 
turns  out  to  loe  the  more  significant  part  of  the  problem.  Our  construction  will  be  tjased  on 
some  variants  of  the  ideas  of  multivariate  B-splines. 


L»t  p  <  m  and  Xat  P  denota  the  projection  of  ■*  onto  tP  ,  i.e. ,  *  •  (P*,y) 

for  *  e  bP  .  If  A  la  a  simplex  in  IfiP  with  varticas  Vg,  ...»  ,  then  the  function 

(3.1)  M(x)  vol|u_p{s  a  A  t  Pt  »  x)  ,  X  a  , 

is  a  piecewise  polynomial  (4),  (8)  of  order  m-p+l  which  is,  up  to  a  constant  (vis., 
vol^CA)  ),  completely  determined  by  the  points  PVg,  PV|||  in  rf*  and  has  its  support 

in  their  convex  hull. 

When  p  -  2,  the  piecewise  polynomial  M  has  any  segment  connecting  any  two  of  the 
projected  vertices  PVq  ,  Pv^  as  a  mesh  line,  and  is  a  polynomial  on  any  connected 

set  not  intersected  by  such  a  mesh  line.  For  many  purposes  (Including  ours),  this  results 
in  too  complicated  a  grid  in  the  plane.  This  cm  be  avoided  if  we  replace  the  simplex  A 
in  (3.1)  by  an  appropriate  m-dlmenslonal  convex  polyhedron  B  and  so  define 

(3.2)  Mg(x)  s-  vol^_p{2  eSsPs-x)  ,  xes^. 

The  function  Mg  so  defined  is  again  piecewise  polynomial  of  total  order  m-p-fl  because 
B  can  be  decomposed  into  disjoint  simpllces  and  therefore  Mg  is  a  sum  of  (slmplicial) 
B-spllnea.  It  has  its  support  in  P(B)  . 

A  more  useful  analytic  description  of  Hg  is  given  by  the  identity  (used  first  in  the 
context  of  slmplicial  B-splines  by  Mlcchelll,  see  [8]) 

(3e3)  /p  ♦<x)  Mg(x)  dx  «  /g  ^(Ps)  dz 

R  ^ 

which  liolds  for  any  test  function  (  ^  8  Cg  ).  T)>e  right  hand  side  defines  Hg  as  a 
distribution  on  .  One  of  the  advantages  of  (3.3)  is  that  it  can  be  used  to  define  Mg 

when  B  is  a  convex  polyhedron  of  dlimnslon  q  <  m  .  The  right  hand  side  is  then 
interpreted  as  a  q-dimensional  surface  Integral.  This  definition  also  malces  sense  when 
P(B)  lias  dimension  <  p  .  In  this  case.  Mg  is  defined  only  as  a  distribution. 

Suppose  then  that  B  is  a  convex  polyhedron  of  dimension  q  4  m  .  As  is  pointed  out 
in  [2],  it  is  easy  to  chedc  the  smoothness  of  Mg  by  using  the  differentiation  formula 

(3.4)  Dp,  Mg  -  -Z  (e-n^>Mg^ 

proved  there.  Here,  the  Bj^  are  the  facets  of  B  ,  i.e.,  the  faces  of  dimension  q-1  , 
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which  make  up  the  boundary  of  B  ,  and  n^^  are  their  respective  normals.  Repeated 

application  of  (3.4)  shows  that  can  be  written  as  a  linear  combination  of 

Mp's  where  each  F  is  a  face  of  B  of  dimension  q  -  a  -  0  .  In  particular,  if  dim  P(F) 

(0,6 ) 

=  p  for  all  faces,  then  D  is  again  a  piecewise  polynomial.  Thus  Mg  is  in 

(,q-d-2  ^  with  d  the  largest  integer  for  which  there  is  a  face  F  of  B  of  dimension 

d  with  P(F)  of  dimension  <  p  .  Note  also  that  the  discontinuities  of  Hg  occur  across 

the  (p-1 ) -dimensional  sets  P(F)  ,  with  F  a  face  of  B  . 

We  now  turn  to  the  construction  of  spline  approximants  on  the  triangular  partition  A 

of  .  Let  e,  i»  (4.  be  the  unit  coordinate  vectors  in  R®  and  define 

t  1]  3=1 

V,  =-  e,  ,  Vj  :=  ej  , 

fe^  +  e^  ,  if  j  «  Mmod  3)^ 

+  e^  ,  if  j  =  2(mod  3)  \  ,  j  =  3,4, ...,m. 


3  “  3(mod  3)  J 

The  vectors  v^ ,  ...,  v^j^  determine  a  parallelepiped  B  with  vertices  s”  ^i''i  ' 

e  {0,1}  ,  all  i  .  The  translated  parallelepipeds  B^^  :=  ie^+je^  +  B  have  pairwise 
disjoint  interiors  and  fill  out  the  slab  x  (0,1)“"^  ,  hence  form  a  partition  for  that 


Consider  the  functions  M  :«  Mg  and  M^j(x,y)  :=  M(x-i,y-j)  »  Mg^  (x,y)  defined  by 
(3.3)  (or,  equivalently,  by  (3.2))  with  p  s  2. 


2.  (i)  E  M^^  »  1  on  R''  . 

(ii)  Mj^j  S  for  r  :»  m-1  and  P  :■» 


Proof.  The  identity  (i)  follows  from  the  representation  (3.2)  and  the  fact  that 
»  R^  X  [0,1]®”^  ,  with  the  having  pairwise  disjoint  interiors. 

We  have  already  noted  that  M  is  a  piecewise  polynomial  of  total  order  m-1  =  r  .  To 
check  the  smoothness  of  M  ,  note  that  any  face  F  of  B  is  of  the  form 

F  »  {v  ♦  E  a.  V.  :  0  <  a.  <  1 } 

•  r  r  1 

r6I 

with  V  some  vertex  of  B  and  I  j  .  if  p(f)  is  a  segment,  i.e.,  p(F)  has  dimension 
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1,  then  all  18  1  are  congruent  aodulo  3  ,  hence  |l|  <  .  Therefore,  M  has 

aeoothneas  p  i»  ■  -  -  2  «  .  Such  a  segeent  P(F)  eust  lie  on  the  line 

P(V)  ♦  XP(Vj)  ,  -■  <  X  <  ••  , 

with  j  •  Keod  3)  ,  all  18  1,  and  j  8  J  .  This  is  one  of  the  eesh  lines  which  fore  the 

partition  t  ,  Since  M  is  discontinuous  only  across  such  segaents  P(F)  ,  it  follows 

that  M  is  in  S  . (A )  ,  and  therefore  so  are  all  the  M,  ^  .  Ill 
r,P  13  '  '  ' 


2r-4  3 

I«Ma  3.  P  >  — ^ —  ,  then  every  s  8  S^  p(A)  can  be  represented  on  as  a 

linear  coebination  of  the  truncated  potiera  in 

T  :•  <“  {***y‘*.  *^<y-yj)^‘*,  (x-y-Xj^)^**(x-»y)^  :  p,q>0  I  P«q<ri 

p,q  >  P  whenever  they  exponentiate  a  truncated  function}  . 


Proof.  Denote  by  0^^  and  the  left,  respectively  right,  triangle  in  A  which 

make  up  the  square  .  Let  a  8  S  and  let  tQ  be  the  polynoeial  of  order  r  which 


agrees  with  s  on 

o*  ,  with 
00  ' 


Than  s  -  tQ  «  0  on  the  line  y  ■>  x  .  nierefora  s  -  tQ  -  u  on 


u(x,y)  -  Cp^(x-y)^(x+y)‘*  . 

P4q<r 

Sot  t,  :«  tQ  +  I  Cpq(X“y)+*’(x+y)**  •  Then  s  ••  t,  on  0^  .  Wo  now  continue  in  this 
manner,  moving  to  the  right,  using  the  truncated  powers  (x-x^)^'Py‘*  and 
(x-y-Xj^)^*’(x+y)^  to  construct  a  function  t  in  the  span  of  T  which  agrees  with  s  on 

m 

P^Q  ,  i-0,1,2,...  .  On  the  first  column  squares,  this  t  is  just  the 

polynomial  tQ  .  Proceeding  now  up  that  column,  we  further  modify  t  by  the  addition  of  a 
suitable  linear  combination  of  the  truncated  powers  x^Cy-yj^)^**  and  (x+y)P(x-y-yj 
(note  that  (x-y-y^)^**  *  (x-y-Xj^)^'*  for  suitable  )  and  so  obtain  a  now  t  in  the  span 
of  T  which  agrees  with  s  on  the  first  row  and  first  columnn  of  squares,  l.s.,  on 

m 

u  a  u  o  » 
i-o  " 

We  claim  that  now  s  -  t  on  all  of  R^^  .  We  prove  this  by  showing  that  s  -  t  on  , 

l,j  >  0  ,  using  lexicographic  ordering  and  induction. 
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Assume  that  s  >  t  on  for  all  ii  <  1  and  for  all  (ti>v)  with  v  »  i  and 

V  <  j  (as  Is  the  case  for  (i>j)  •  (1>1)  )•  Then  u  s-t  Is  a  piecewise  polynomial 
function  of  total  order  r  and  smoothness  P  which  vanishes  for  0  <  x  <  and  for  0  < 
y  <  .  If  now  u  were  nonzero  somewhere  In  ,  then,  on  restricting  u  to  some 

lino  x+y  »  c  ,  with  Xj^+y^  <  c  <  *i+i+yj  »  wo  would  obtain  a  univariate  piecewise 
polynomial  function  v  of  order  r  and  smoothness  P  which  vanishes  outside  some 
interval  [a,b]  and  has  just  three  knot  locations,  viz.  the  points  a  ,  (a-fb)/2  ,  and  b  . 
Further,  v  would  be  nonzero  somewhere  in  ta,b]  .  This  would  imply  that  the  sum  of  the 
multiplicities  of  the  knots  a  ,  (a-fb)/2  ,  and  b  is  at  least  r+1  ,  which  would  imply 
that  at  least  one  knot  has  multiplicity  <  ,  therefore  P  <  r-1  -  ™_| 

,  a  contradiction.  This  advances  the  induction  hypothesis  and  so  finishes  the  proof.  | | | 

Slieozem  3.  us  (d  )  is  dense  in  if  and  only  if  P  <  . 

n  r,p  n  -  w  '  —  * -  3 

Froof.  By  Lemma  2,  S  .(&)  contains  the  local  and  positive  partition  of  unity 
r  ,P 

(M,j)  as  long  as  P  <  .  This  implies  the  density  of  U  S  „(d  )  by  the  following 

ij  3  nr^pn 

standard  argument;  The  simple  approximation  map 

Tf  Z  f(Xj,yj)Mj^j 

carries  f  into  satisfies 

I  (f-Tf  )(x,y)I  =•  I  Z  (f(x,y)-f(x  ,y.  ))m  (x,y)| 

i,j  ^  ’ 

<  max  |f(x,y>  -  f(x  ,y  )| 

M.^(x,y)^0  J 

Consequently,  If  -  Tfl  <  u(f;  diam  supp  N)  .  Now  scale. 

2r-4 

For  the  converse,  assume  that  P  >  — —  .  Consider  the  difference  operator 
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A.  t»  .  Since  Af  and  Af.  .  .  coemute,  %#e  have  A.  t  •  0  for  all  t  In  the 

h  h  (hfh)  h  (h#n;  h 

2 

span  of  defined  in  tjerama  3  .  Therefore,  by  Lemma  3,  A^s  •  0  on  for 

every  s  6  S  _ (A )  . 

f  »P 


“,_2 


Suppose  now  that  f  e  C^CR  )  and  that  If  -  s^l  0  (as  n  ♦  for  some  s^  6 

2 

S  (A  )  ,  na1,2,...  .  Assume  without  loss  of  generality  that  f  has  Its  support  In  H.. 
r,P  n  ^ 

Since  A^  is  obtained  from  A  by  scaling,  it  follows  that  »  0  on  R^^  ,  therefore 

A^f  •  0  .  Dividing  by  h^*^  and  taking  the  limit  as  h  0  shows  that 


(®e  “e  “e  ie  °  ' 

“l  ®2  *1^2 

•  2  III 

idilch  shows  that  f  la  not  an  arbitrary  function  In  C.(R  )  .  |  |  | 


Next  we  consider  the  approximation  from  S  .(A  )  to  smooth  functions.  For  this,  w 

r  ,p  n 

need  to  consider  the  B-spllnes  associated  with  faces  of  B  .  Any  X  ^  B  Is  associated  with 
a  face,  viz.  the  face 

F  F.  {  r  O.v.  :  0  <  o  <  1}  . 

lei  ^ 

This  Is  a  face  of  dimension  | I |  ,  or,  a  | I | -face,  for  short. 

Denote  the  corresponding  B-spline  by  .  We  are  particularly  interested  in 

faces  for  which 


2  c  I  and  I  >  2 


For  such  a  face,  set 


0  «-  {  I  <»,v  s  0  <  o  <  1}  . 

101 

Then,  for  any  test  function  ♦  , 


/  ,  -  /p  ♦op  , 


using  the  facts  that  ♦(Pz)  -  ♦(Pz^)  in  case  z  6  z^  +  Q  ,  and  that  the  (m-|ll)- 
dlmenslonal  volume  of  Q  Is  1  .  This  shows  that,  for  such  a  face, 

(3.5)  M^(x,y)  -  vol|j_2{z  9  F+Q  :  Pz  -  (x,y)>  -  ''°^|l|-2^“o  '  ***0  ”  • 


Now  define 


”i1  '*  "f  ' 
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Lat  now  8  >  3  ,  assuae  th«  clala  provad  for  all  | I |  <  a  ,  and  consldar  soaa  I 
with  |l|  -  8  . 

If  d  -  |l|  -  1  ,  then  h  >  1  and  tha  clala  radueaa  onca  again  to  (3.5).  Otherwlae. 
d  <  |l|  -  1  .  In  thl8  caaa,  conaldar  |i  8  I  \3  >od  aet  u  .  He  calculate 

for  some  0+8  <  h  .  By  (3.4), 

(3.9)  Dp^  -  -J:(u.n^) 

where  the  are  the  ( ll (-1 )-face8  of  Fj  and  the  n^  ate  the  corresponding  outward 

noreals.  this  lepllea  that  (u*n^)  -  0  for  any  face  F^  parallel  to  u  .  There  are  only 
two  faces  «(hlch  are  not  parallel  to  u  ,  vis.  the  faces  Fj  and  Fj+u  ,  with 

J  I*  I  \  (u)  . 


Since  their  noreals  sue  to  sero  and  |(u*n^)|  •  1  ,  (3.8)  becoeea  sleply 


This  lepllea  that 
(3.10) 


-  +  H? 


“Pu^«8  -  -li 

J 


With  I-  X  y  -  C  n  and  (C#n)  *-  <x,y)  -  Pu  •  Kora  explicitly. 

If  Pu  -  /  Sj  \  then  fa0(i»j)  ■<  Bi  ^  >  • 

le^+ejJ  al®"’j®+  8l“j®"j 

Since  hj  >  h  -  1  and  |j|  <  |l|  and  f^^  Is  of  coordinate  degree  <  (a, 8)  ,  we  nay  apply 
the  Induction  hypothesis  to  conclude  that  •  polynoelali  eore  precisely,  that 


°Pu^<.8‘*'^’ 


0-1  8 
ox  y 

a  O  8-1 

8x  y 

0-1  8  -  o  8-1 
ox  y  +  8x  y 


+  lower  order  teres.  If  Pu 


Integrating  bacic  up,  wa  find  that 

a  B 

♦o8^*'^*  -  X  y  +  lower  order  teres  ♦ 

with  g  a  univariate  function  and  s  8  perpendicular  to  Pu  . 

But  now,  since  |l|  >3  and  d  <  |l|  -  1  ,  It  Is  possible  to  choose  u  In  different 
ways,  say  u  •  v  and  u  •  w  ,  so  that  Pv  and  Pw  are  linearly  Independent.  This  allows 


the  conclusion  that  actually  (3.8)  holds  and  so  advances  the  Induction  hypothesis. 
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This  finishes  the  proof  of  the  Clale*  Tetlclng  now,  in  particular,  I  •  g  ,  we  find 
that  d  "  '  hence  hj  •  m  -  “  L'^— i  “  **  '  “*  conclude  that  (3.8)  holds  for 

all  a-fd  <  )c  .  Repeated  application  of  this  fact  gives  us,  for  each  O'fB  <  )c  ,  a 
polynoeial  of  coordinate  degree  (a, 8)  with  leading  coefficient  1/(alBl)  such 

a  $ 

E  Q,g<i,j)  M^^(x,y)  -  .  Ill 

*•  J 


We  conclude  that,  for  any  r  ,  e  ,  and  any  O40  <  k.  , 

.a,..  -tP 


(x-r)  (y-8)  „  ^  . 

- Jril -  ■  /  Vr,J-s‘*-''y-*> 

J 


E  Q^g(i-r,J-#)  M,^(x,y) 


i.J 


ij 


Consequently,  for  any  p  6  , 


01  3 

(3.11a)  p(x,y)  -  E  (D‘“'®’p)(r,s)  -  E  q„(i,j)  M  (x,y) 

o+fl<k  i,j  ” 

with 

(3.11b)  ^  Q„a(*-r»y”«)(D*’*'®'p)(r,s)  . 

o+B<)c  ’’ 


The  next  lesuna  implies  that  is,  in  fact.  Independent  of  r  and  s  . 


S.  The  map  ^  :  ql — ►  E  q(i,j)  is  one-one. 


i.j 


Proof.  Suppose  q  6  P)^\{0}.  Then,  for  any  r  >  0  ,  there  exists  c  9  so  that  q 
is  of  one  strict  sign  on  the  ball  B^(c)  of  radius  r  and  center  c  .  Choose  r  = 
diam  supp(H)  .  Since 

(i,j)  8  supp(Mj^j)  -  (i,j)  +  supp(M)  , 

it  follows  that  {q(i , j )M^^ (c )  :  M^^(c)  *  0)  are  all  of  one  strict  sign,  hence 

E  q(i,j)M^^  *  0  .  Ill 


l.j 


He  conclude  that  q 


in  (3.11)  does  not  depend  on  r  ,  s  .  In  particular 


-  qt.(i.J)  -  I  Q_6(0,0)(d‘*'®’p)(1,J)  . 

a+S<k 

Oorollaty.  For  any  p  fl  ,  p  •  E  ^(p(*+l» •+! ) )  M. .  with 

i.J  ^ 

'■  ^  (D‘“'®’f)(0,0)  • 

a+e<lc 

Ne  use  the  linear  functional  1  to  define  an  approxlaatlon  fro«  p(A)  In  the  now 
standard  quasl-lnterpolant  fashion.  First,  we  modify  X  so  as  to  make  It  applicable  to 
any  f  8  .  For  this,  let  v  be  a  bounded  extension  of  X  from  to  C(o^*)  . 

Th^n  u  can  be  baken  to  be  a  bounded  linear  functional  on  all  of  CqCK  )  •  In  this  way#  we 
obtain  a  linear  wap 

Lf  t-  E  Jl(fC+l,*.+j))  M 

,  i'J 

on  CnO^)  to  S  (A)  which  la  local,  reproduces  P^  ,  and  Is  bounded  by  llii  since 
o  r,P  * 

^Ij^lj  ~  ^  Implies  that 

(3.12)  If  -  bfl(A>  <  dill  +  1)  dlstjjj^j(f,  P^) 

with  N(A)  U  {supp(Mj^^)i  aupp(Mj,^)nA  *  (»}  .  Scaling  by  1/n  to  get  to  the  partition 

A  gives  the  map  L„  to  S  ^(A_)  for  tritlch  we  have  the  following 
n  n  r,p  n 

thmocea  4.  Set  m  r-1,  P  «-  and  X  «-  L^f^J  -  P+2  •  il  «  « 

then  L_f  0  S_  .(A_)  and 
-  n  r,P  n  - 

(3.13)  If  -  L^fl  <  constj^  '  "“’i*,... 

with  Wj^  the  k-th  order  modulus  of  smoothness  of  f  .  In  particular.  It  f  8  C^(*^)  ,  then 

(3.14)  If  -  Ljjfl  <  constj^  |f|j^  n“’‘  . 


Proof.  l<et  A 


By  (3.12), 


o^^(n)  and,  correspondingly 

N(A)  ;«  U{supp(M  ):  supp(M^_)nA  *  ft)  . 
pqn  pqn 


If  -  L  fl(A)  < 
n 


(const  +1)  dlst„,,. 
r  N(A) 


(3.15) 


(f,  P^) 


vihll*,  c.g*,  froa  [5] 


With  c  m  W|^(f,  1/n)  or  c  ■  ainc*  dlaa  MCA)  ■  0(1/n)  •  ||| 

4«  Ceaoladiag  rwaarka*  Tha  eonatructlor.  of  aaooth  auXtivarlata  apllna  Intarpolanta 
davalopad  In  Saetlon  3  can  ba  axtandad  to  aora  ganaral  parti tlona  (aaa  alao  (3])«  Suppoae 
that  n  la  a  partition  of  which  can  ba  obtalnad  aa  followa.  Starting  with  a  regular 
partition  11^  ,  aaaoclated  with  the  two  independent  dlreotlona  d^  and  ^2  •  *** 
llnea  through  all  the  vertleea  of  0^  In  the  dlreotlona  of  the  veetora  dj,  dp  , 

with  these  dj^'s  nonsero  vertices  of  11^  other  than  d^  or  dj  >  and  palrwlae  linearly 
Independent. 

Thus  I  la  aaaoclated  with  the  directions  a^  and  ej  >  A  is  associated  with  the 
directions  e,,  e^  and  ,  while  the  four  directions  e^,  Sj. 

aaaoclated  with  a  partition  In  %ihlch  both  diagonals  ara  dratm  Into  every  square. 

Given  such  a  partition  n  and  n  >  p  ,  define  vectors  v^,  v_  In  S*  by  v^  :» 

d^  /  Vj  >>  dj  /  and  ,  for  1  >  2  ,  v^  dj  +  e^  ,  with  j  6  p  and  j  ■  Knod  p)  . 

Set  B  >•  {£*  o^v^  <  0<a^<l}  .  then  the  corresponding  B-spllne  Mg  given  by  (3.3)  la  of 
order  r  i"  b-1  and  anoothneas  P  i"  r-1  -  “f*  the  vertleea  of 

,  than  the  translated  splines  given  by  >>  Mg(x-x^,y-y^)  ,  all  1  ,  fom  a 

local  partition  of  unity.  Using  arguaents  like  those  for  Lema  4  and  Theores  4  ,  they  can 

ba  used  to  construct  a  local  and  bounded  quasllnterpolant  on  Cg(S^)  Into  p^^n^  which 

k  2 

reproduces  polynowlals  of  total  order  k  >•  9*2  ,  hence  approxieates  )-functions  to 

within  0(n“’')  . 

The  particular  choice  dj  ><•  (cos  ^  ,  aln  ,  ja1,2,3,  results  In  a  partition 

of  Into  equilateral  triangles.  He  have  recently  learned  that  P.  Frederlckson  161,  [71 

has  studied  spline  approxlsatlon  on  this  partition  and  has  shown  the  existence  of  spline 

Intarpolanta  frow  S  _(n  )  In  the  case  p  ■  ^r^d  ^  _  2(Bod  3).  This  is 

r,p  n  3 

accospllshed  by  using  a  partition  of  unity  given  as  translates  of  a  fixed  finite  support 
spline  obtained  by  a  certain  convolution. 
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